Green's function for second order parabolic 
systems with Neumann boundary condition* 
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Abstract 

We study the Neumann Green's function for second order parabolic 
systems in divergence form with time-dependent measurable coefficients 
in a cylindrical domain Q = Q x ( — 00, cc), where f2 cz R" is an open con- 
nected set such that a multiplicative Sobolev embedding inequality holds 
there. Such a domain includes, for example, a bounded Sobolev extension 
domain, a special Lipschitz domain, and an unbounded domain with com- 
pact Lipschitz boundary. We construct the Neumann Green's function in 
Q under the assumption that weak solutions of the systems satisfy an inte- 
rior Holder continuity estimate. We also establish global Gaussian bounds 
for Neumann Green's function under an additional assumption that weak 
solutions with zero Neumann data satisfy a local boundedness estimate. 
In the scalar case, such a local boundedness estimate is a consequence of 
De Giorgi-Moser-Nash theory holds for equations with bounded measur- 
able coefficients in Sobolev extension domains, while in the vectorial case, 
one may need to impose further regularity assumptions on the coefficients 
of the system as well as on the domain to obtain such an estimate. We 
present a unified approach valid for both the scalar and vectorial cases 
and discuss some applications of our results including the construction 
of Neumann functions for second order elliptic systems with measurable 
coefficients in two dimensional domains. 
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1 Introduction 

In this article, we are concerned with Green's function for parabolic systems with 
Neumann boundary condition, which hereafter we shall refer to as Neumann 
Green's function. It is often called Neumann heat kernel if the parabolic system 
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has time-independent coefficients. In a cylindrical domain x (— cx),oo), where 
the base is an open connected set in M", we consider parabolic differential 
operator of the form 

Here, we use the summation convention over repeated indices 1 ^ a, /3 ^ n, it 
is a column vector with components u^, . . . , u^, A°'^{x, t) arc N x N matrices 
whose elements (a;, t) are bounded measurable functions defined on the entire 
space R"+^ such that for arbitrary N x n matrices ^, rj with elements 77^, 
we have 

\\e < < ^"'i^ii^i' (1-2) 

where A is positive constant and = Ca^a- We emphasize that we allow 
coefficients to be non-symmetric and time-dependent. By a Neumann Green's 
function for the system (|l.ip in il x (—00, 00), we mean a.n N x N matrix valued 
function N(x, t, y, s) that is a solution of the initial-boundary value problem 

J^x,tN{x,t,y,s) = Q in f7x(s,cx)), 

^N{x,t,y,s) = on cQ. x (s,cx)), 
N{x,t,y,s) = Sy{x)I on x {t = s}, 

where o/cj/ is the natural co-normal derivative, 6y{-) is Dirac delta function 
concentrated at y, / is the N x N identity matrix, and the equation as well as 
the boundary condition should be interpreted in some weak sense. Equivalently, 
it can be defined as a function that weakly satisfies 

j!fj:^tN{x,t,y,s) = Sy{x)Ss(t)I in f}x(-a3,co), 
■^N{x,t,y,s) = on c?fi x (-cx), <»), 

N{x,t,y,s) = on fix(— cx),s). 

More precise definition of Neumann Green's function is given in Section [2.41 In 
the case when the coefficients are time- independent, K{x, y, t) ■■= N{x, t, y, 0) 
is usually called a Neumann heat kernel. 

We prove that if the base ft is such that the multiplicative embedding in- 
equality p.ip holds there and if weak solutions of the system (jl.ip satisfy an 
interior Holder continuity estimate (|3.2p . then the Neumann Green's function ex- 
ists and satisfies a natural growth estimate near the pole; sec Thcorcm l3.9l The 
multiplicative embedding inequality p.ip holds if, for example, f2 is a Sobolev 
extension domain with finite Lebesgue measure, a special Lipschitz domain, and 
an unbounded domain with compact Lipschitz boundary, etc. Also, the Holder 
continuity estimate p.2[) holds, for example, in the scalar case (i.e., N = 1), in 
the case when n = 2 and the coefficients are time- independent, and in the case 
when the coefficients belong to the VMO class; see Section |4l 

We are also interested in the following global Gaussian estimate for the 
Neumann Green's function: For any T > 0, there exists a positive constant C 
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such that for all t,s satisfying < t — s < T and x, ?/ e fi, we have 

m-,t,y,s)\ < (^T^expj-^^l , (1.3) 

where k is a positive constant independent of T. We present how one can derive 
the global estimate like (|1.3p using a local boundedness estimate (|3.18p for weak 
solutions of the system with zero Neumann data; sec Theorem 13.211 It is well 
known that the local boundedness estimate p.l8|) is available, for example, when 
A'^ = 1 and is a bounded extension domain or when n = 2, the coefficients are 
time- independent, and 17 is a bounded Lipschitz domain, etc.; see Sectional In 
fact, we show that the local boundedness estimate is a necessary and sufficient 
condition for the system to have a global Gaussian estimate for the Neumann 
Green's function; see Theorem 13.241 

The Green's function for a parabolic equation (i.e., = 1) with real measur- 
able coefficients in the free space was first studied by Nash [53] and its two-sided 
Gaussian bounds were later obtained by Aronson [2] . There is a vast literature 
on Green's function for parabolic equations satisfying Dirichlet or Neumann 
conditions. In the case when the coefficients are time- independent, they are 
called Dirichlet or Neumann heat kernel and have been studied by many au- 
thors; see Davies [12] and references therein, [3S1 |3I] for a treatment of heat 
kernels on non- Euclidean setting, and also a very recent article [18]. It is well 
known that Aronson's bounds are no longer available for a parabolic equation 
with complex valued coefficient when n > 3. Related to Kato square root 
problem, Auscher [3] obtained an upper Gaussian bound of the heat kernel for 
an elliptic operator whose coefficients are complex L"-perturbation of real co- 
efficients; see [H [SJ |B] and also [TJ [TTl [5S] for related results. The Dirichlet 
Green's function for a parabolic system with time-dependent, measurable co- 
efficients was studied in papers by Cho et al. [319], where Dirichlet Green's 
function is constructed in any cylindrical domain under an assumption that the 
weak solutions of the system has an interior Holder continuity estimate that 
is equivalent to ours. In many respects, this article can be considered as nat- 
ural follow-up of their papers. However, we must point out it requires some 
effort and caution to handle Neumann condition and there are certainly some 
subtleties involved here. The novelty of our paper lies in the following points. 
First, we allow time-dependent, non-symmetric coefficients and do not specifi- 
cally assume any regularity on them. Secondly, we treat parabolic equations and 
systems simultaneously and do not rely on methods that are specific to single 
equations case, such as the maximum principle or Harnack's inequality. Finally, 
our method separates existence of Neumann Green's function from its global 
estimates, which enables us to construct Neumann Green's functions in a large 
class of domains, where global Gaussian estimates for Neumann functions may 
not be available anymore. In this paper, we only deal with so called strongly 
elliptic systems as it is described in (|1.2p , but our method can be as well applied 
to systems satisfying a weaker condition; see [29] for the Neumann heat kernel 
for the elliptic system of linear elasticity. For the heat kernel for second-order 
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elliptic operators in divergence form satisfying Robin-type boundary conditions, 
we bring attention to a very recent article |16j . 

The organization of the paper is as follows. In Section [51 wc introduce some 
notation and definitions including the precise definition of Neumann Green's 
function. In Section |31 we state our three main theorems, which we briefly 
described above. In Sectional we give examples for which our main results apply, 
and also in Theorem 14. 5 [ we construct the Neumann function for second order 
elliptic systems with measurable coefficients in an extension domain $7 c 
and obtain logarithmic bound for it. We give the proofs for our main results 
and Theorem 14.51 in Section [SJ 

2 Preliminaries 

2.1 Basic notation and definitions 

Wc use X = {x,t) to denote a point in ]R"+^; a; = (xi, . . . ,x„) will always be 
a point in M". We also write Y = (y, s), Xq = {xo,to), etc. We define the 
parabolic distance between the points X = {x, t) and Y = (y, s) in R"+^ as 

\X -Y\,^ := max(|a; - y\ , s/Jt^^l) , 

where | • | denotes the usual Euclidean norm. Wc write \X\g^> = \X — 0\,^. For 
an open set Q c we denote 

dx = dist{X, dpQ) = M{\X -Y\^: Y e dpQ}; inf = cx), 

where dpQ denotes the usual parabolic boundary of Q. We use the following 
notions for basic cylinders in ]R"+^: 



where Br{x) is the usual Euclidean ball of radius r centered at x e M". We use 
the notation 



We denote a a b = min(a, b) and a v b = max(a, b). 
2.2 Function spaces 

For spaces of functions defined on a domain ft c K", we use the same notation 
as used in Gilbarg and Trudinger [T7], while for those defined on a domain 
Q d K^+^j we borrow notation mainly from Ladyzhenskaya et al. [22]. To 
avoid confusion, spaces of functions defined on Q c ]R"+^ will be always written 
in script letters throughout the article. For q ^ 1, we let JifqlQ) denote the 



Qr{X) 



Br{x) X (t - r'^,t), 
Br{x) X {t,t + r'^), 
Br{x) X {t-r'^,t + r^) 
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classical Banach space consisting of measurable functions on Q that are 5th- 
integrable. ^q^r{Q) is the Banach space consisting of all measurable functions 
on Q with a finite norm 

where 9 5= 1 and r ^ 1. ^9,9 (Q) will be denoted by ^q{Q). By ^^'^/^(Q) we 
denote the set of all bounded measurable functions w on Q for which |u|p,^/2;Q 
is finite, where we define the parabolic Holder norm as follows: 

I"UW2;Q = Ma'W2;Q + l'"lo;Q 

== sup 'ly^^y + snpKX)|, /.6(0,1]. 

We write u e ^"(Q) (rcsp. 'i^^{Q)) if w is an infinitely differentiable function on 
with a compact support in Q (rcsp. Q). We write DiU = Dx^u = du/dxi 
(i = 1, . . . , n) and Ut = du/dt. We also write Du = D^u for the column vector 
{Diu, . . . , Dnu)^ . We write Q{to) for the set of all points {x, to) in Q and I{Q) 
for the set of all t such that Q{t) is nonempty. We denote 

III'^IIIq ~ / \Dxu\'^ dx dt + CSS Slip / \u{x,t)\'^ dx. 

The space #^^'°((5) denotes the Banach space consisting of functions u e ^q{Q) 
with weak derivatives DiU e ^q(Q) (i = 1, . . . , n) with the norm 

ll"ll<'°(Q) = ll"ll^,(Q) + \\Dxu\\j^^(^Q) 
and by i^^'^iQ) the Banach space with the norm 

ll"ll<''(Q) = ll"ll^,(Q) + \\^M\x^{Q) + l|wt||.S-,(Q)- 

In the case when Q has a finite height (i.e., Q a K" x (— T, T) for some T < <»), 
we define y-iiff) as the Banach space consisting of all elements of #^^'°(Q) 
having a finite norm ||'u||y2(Q) •= III'^IIIq and the space 'f'^^'^iQ) obtained by 
completing the set #^^'^(Q) in the norm of "^(Q). When Q has an infinite 
height, we say that u e r2(Q) (resp. r2^'°(Q)) if u e f2{QT) (resp. ■^■^'^{Qt)) 
for all T > 0, where Qt = Q n {\t\ < T}, and |||m1||q < 00. Note that this 
definition allows that 1 G ■>^^'°(f2 x (—<»,(»)) when |n| < 00. Finally, we write 
u 6 ^q^iociQ) if u 6 J^q{Q') for all Q' (szQ and similarly define '^qfociQ)^ etc. 

2.3 Weak solutions 

For / 6 -Sf2,i(Q)^ and e ^2{Q)'^ (a = 1, . . . ,n), we say that m is a weak 
solution of = f Dag"" in Q if tt 6 "^(Q)^ and satisfies the identity 

M /" i^(M,0)= / /•</>+ /" (2.1) 

I JQ JQ 
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for all e "S^iQY , where we set 

The adjoint operator ^* of ^ is given by 



Note that the coefficients of p.2p also satisfy the ellipticity and boundedncss 
condition ()1.2p with the same constant A. We say that it is a weak solution of 
= / - D^g" in Q if M 6 r2(g)^ and satisfies, for all e '^;r^'(Q)^, the 
identity 

= /" / / (2.3) 

Let fl be an open set in R" and E be a relatively open subset of c?J7. Below we 
denote Q = fl x {a,b) and S' = Ex (a, 6), where — cx) ^ a < 6 ^ <». We say that 
u is a weak solution of 

= f Dag" in Q, du/dv = g°'i'a on S 

if M 6 r2(g)^ and satisfies the identity ^ for aU e "^f (Q u S)^ . Similarly, 
we say that it is a weak solution of 

= / - Dag°' in g, (5M/(5i/* = g°'va on 

if M 6 fziQ)^ and satisfies the identity (gS]) for all e '^^{Q u 5)^. For 
■0 6 and a 7^ — cX), we say that m is a weak solution in ^^^'^(Q)^ of the 

problem 

' ^u = f- D^g" in Q 
du/dv = g^'i^a on dpQ 
u{-, a) = ip on ri 

if M 6 y^'^iQ)^ and satisfies the identity 

JQ JQ JQ JQ Jn 

for all cf) e "^^^(Q)^ that is equal to zero for t = b. For b ^ 00, wc say that u is 
a weak solution in 7^^'°((5)^ of the problem 

= f D^g"^ in Q 
du/du* = g^Va on dpQ 
u{-, b) = tp on i7. 

if M 6 'f'2'^{Q)^ and satisfies the identity 

/ M-0,+ /" ^(0,m)- / f ct>- [ g"-Dact>= [ ^■ct>{;b) 

JQ JQ JQ JQ Jn 

for all 4> e ^S/f^iff)^ that is equal to zero for t = a. 
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2.4 Neumann Green's function 

Let Q ■■= il X (—<»,(»), where f2 is a domain in R". We say that an N x N 
matrix valued function N{X, Y) = N(x, t, y, s), with measurable entries 

Q X [-«,«)], 

is Neumann Green's function of ^ in Q if it satisfies the following properties: 

a) For aU Y e Q, wc have N{-,Y) e W^f^^iQ)^" and N{;Y) e mOxQAY))^" 
for any r > 0. 

b) For all r 6 Q, we have ^N{-,Y) = Syl in Q, ^N{-,Y) = on dpQ, in 
the sense that for any (f) e "^^^^(2)^, we have 

1^ {-N,k{;Ym + al^D^N,k(.;Y)Dacj^^) = /(F). (2.4) 

c) For any / = (/I, . . . , f^)^ e '^;^{Q)^, the function u given by 

:= / N{Y,Xff{Y)dY 

Jq 

is a weak solution in r/'"(Q)^ of = f in Q, dujdv* = on 5pQ. 

We note that part c) of the above definition gives the uniqueness of a Neumann 
Green's function. 



3 Main results 

We make the following assumptions to construct the Neumann Green's function 
in Q := X (—00, {»). 

Al. We assume that C^(f2) is dense in Vl^^'^(O) and there exists a constant 7 
such that for all u e W^'^{il), we have 

< l\\Du\\%^l^^'^\\u\\%%'\ (3.1) 



where we use the notation 



w^'^{n) if \n\ = oD 

{ueW^-^{n): J^u = 0} if < cx). 



A2. There exist fiQ e (0, 1], i?c e (0, oo], and Bq > such that if it is a weak 
solution of ifw = (resp. = 0) in Q = Q]^{X) (rcsp. Q = Qr{X)), 

where X e Q and < R < dx ■= dx ^ Rc, then u is Holder continuous in 
= Qfl,/2(^) (resp. = Qfl/2("''^)) ^^^'^ estimate 



1/2 

(3.2) 
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Remark 3.3. In the case when n > 3, the mequahty 

||u||i2n/(„-2)(0) C||L'u||i2(Q) (3.4) 

imphes (via Holder's inequality) the inequahty p.ip and also the inequality 

< C\\Du\\%%;;'^\\u\\f-^^^\ (3.5) 

If f2 = R" . then p.4p and p.5p arc the Sobolev inequality and the Nash inequal- 
ity respectively. In fact, they arc equivalent in that case; see the monograph 
by Saloff-Coste |28] . From this point of view, one may think Al roughly means 
that behaves like M" with respect to W^''^ functions. 

Remark 3.6. The inequality p.ip implies the following: If m G '^(fi x (a, 6)) is 
such that u{-,t) e W^''^{ft) for a.c. t e (a, 6), where — cx) ^ a < b ^ co, then we 
have (see [11 pp. 74-75]) 

ll"ll^2(„ + 2)/„(nx(a,h)) 7|||u|||nx(a,fc)- (3.7) 

Another important consequence of (|3.H) is that for u e 1^(il x {a,b)) with 
6 — a < 00, we have 

ll"lb2„.+2,/„(ax(a,6)) < (27 +{b- a)"/2|f7|-i) ~ |||u|||nx(a,;,). (3.8) 
We refer to [H Eq. (3.8), p. 77] for the proof of H^ . 

Theorem 3.9. Assume the conditions Al and A2. Then there exists a unique 
Neumann Green's function N{X,Y) = N{x,t,y, s) of ^ in Q ■■= ilx (— cx),cx)). 
It is continuous in {(X,Y) e Q x Q: X ¥= Y} and satisfies N(x,t,y, s) = for 
t < s. For any f = (/\ . . . , /^)^ e "^^^(2)^; ^'le function u given by 

u{X) := / N{X,Y)f{Y)dY (3.10) 
JQ 

is a weak solution in '^^''^(Q)^ of 

^u = f in Q, du/8iy = on 8pQ. (3.11) 
Moreover, for all ■0 = {ip^ , . . . , ip'^)'^ e L'^{fl)^ , the function given by 

u{x,t)= [ N{x,t,y,s)iP{y)dy (3.12) 
Jn 

is a unique weak solution in "^^'"^(rj x (s,co))^ of the problem 

J^u = in fl X {s, co) 
< du/di^ = ondnx{s,rx,) (3.13) 
it(-, s) = on ^l. 
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Furthermore, for X,Y e Q satisfying < \X — Y\^0g < ^dy, we have 

\N{X,Y)\^C\X-Y\-^", 



(3.14) 



and for X, X',Y e Q satisfying 2\X -X'\^ <\X - Y\,^ < ^dy, 

\N{X,Y)- N{X',Y)\ ^ C|X-X'|^|X-y|^"~^", (3.15) 



where C = C{n, N, A, /io, Bq). 

Remark 3.16. In the proof of Theorem 13.91 we will sec that for any Y e Q and 
< R < dy, wc have 



|||iV(-,F)||ls\Q,(y)<Ci?-"/2. 

n) ||A^(-,r)||^^^„^^^^„(2\Q^(^)) 

Hi) \{X e Q: \N{X,Y)\> t}\ ^ C-yT-"^, Mrxdy'^. 

iv) |{Xe Q: |D^Ar(X,y)| > t}| Cr"-^, Vr > dy^"^^^ 
v) ||Ar(,r)||^^(Q„(^))<C,,,i?-^+("+2)/p^ Vpe[l,H±2). 

vz) p,Ar(.,y)||^^(Q^(^)) < Vpe [1,^^). 

Here, we denote 7V(-, Y) = N{-,Y) - jkHit - s)/, where i7(t) = l[o,x.)(*) and 



we use the convention l/oo = 0. 

Remark 3.17. If fl a R" is such that there exists a bounded trace operator 
from W^^'^in) to L^{^), then it can be shown that for / e '^[^'{Q x [a, b])^ and 
g e '^^°°(5ri X [a, 6])^, where — oo < a < b < co^ u defined by 



The following "local boundedness" assumption is used to obtain global Gaus- 
sian estimates for the Neumann Green's function. 

A3. There exist Rm e (0, go], i?i > such that if it is a weak solution of 



= in Q = Qjf{X) n Q, du/dv = on 5 = Q^{X) n dpQ, 
(resp. = in Q = Q\{X) n Q, dujov* = on S" = Qj(X) n o^Q, ) 




is a unique weak solution in y^ ^i^ x (a, b))^ of the problem 

' ^u = S in f7 X (a, 6) 
< du/dv = g on pfi x (a, 6) 
«(•,«) =0 on ri. 
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Remark 3.19. Note that A3 implies that there exists > such that for all 
X e il. and < i? < Rm, we have 

OR'' \n n Br{x)\, (3.20) 

which follows from (|3.18p by taking u = 1 and 6 = \/B\. In particular, we 
should have Rm < co in the case when \U\ < co. In fact, we may take Rm = 
in that case without loss of generality. 

Theorem 3.21. Assume that the conditions A1-A3. Then, for t > s, we have 
the Gaussian hound for the Neumann Green's function 

where C = C(n, N, A, Bi) and k = k(A) > 0. In particular, the estimate (|3.22p 
implies that for < t — s ^ T, we have 

\N{x,t,y,s)\^—^cxp(^^^^^^], (3.23) 



{t - s)"/2 t-s 
where C = C{n,N,\Bi,T). 

Finally, the next theorem says that the converse of Theorem 13.211 is also 
true, and thus that the condition A3 is equivalent to a global bound p.22p for 
the Neumann Green's function. 

Theorem 3.24. Assume that the conditions Al and A2. Suppose there exist 
constants B2 and k > such that the Neumann Green's function has the hound 

mx.t,y,s)\ ^ — ^^^exp ( , (3.25) 



{{t-s)^Ri,ri' \ t-s , 

for all t < s and x,y e fl. Then the condition A3 is satisfied with the same Rm 
and Bi depending on n, N, A, k, and i?2. 

4 Applications 

4.1 Examples 
4.1.1 Examples for Al 

First, we give examples of domains satisfying the condition Al. 

(1) We say that 17 is an extension domain (for W^''^ functions) if there exists a 
linear operator E: W^''^{n) W'^^'^iR") such that 

\\Eu\\l2(^^^ (^111411^2(0), |l£'u|jvi/i.2(R,>) sg C||u||h/i.2(o). (4.1) 



10 



We recall that the multiplicative embedding inequality p.ip holds for all 
u G W^-'^(W^) with n^l; see [22j Theorem 2.2, p. 62]. If 1] cz M" (71 5= 1) 
is an extension domain satisfying < co, then for u G VF^'^(f2), we have 

where wc used the Poincarc's inequality in the last step, namely. 

Therefore, we have Al in this case when c R" (n ^ 1) is an exten- 
sion domain in with finite Lebesgue measure. Extension domains include 
locally uniform domains considered by Jones |20] : the extension operator 
that originally appeared in [20] may not satisfy the first inequality in (|4.ip . 
but Rogers [57] constructed an extension operator that satisfies both in- 
equalities in (|4.ip . 

(2) Let fi be a special Lipschitz domain in R" {n 5= 2), i.e., a domain above 
a Lipschitz function ip: R"^^ R. Then, there exists a linear extension 
operator E: W^^^{n) W^-^{R") such that (see [H]) 

\\Eu\\l^R'-) ^ C||u||l2(ji), \\D{Eu)\\l2(^^„) ^ C||Du||L2(a). 
Therefore, similar to the above, for any u e W^''^{n), we have 

More generally, if f2 = U^^j^ ^i, where fli is special Lipschitz domain (ro- 
tating the coordinates axes if necessary), then it also satisfies Al. Indeed, 

h|L.<„..v„(ao < C\\Du\\%%^^f\\u\\%^-^^f ^ C\\Du\\li^;^^'^\\u\\%%'\ 
and thus we have 



(3) Suppose that fl is an unbounded domain in R" {n 5= 3) with compact 
boundary and that fl has the cone property. Then, for all u g W^'^{ft), we 
have (see [7] Theorem 1]) 

||u||i2„/("-2)(o) C\\Du\\L^n), 
which together with Holder's inequality imply Al; see Remark 13.31 
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4.1.2 Examples for A2 

Next, we give examples when the condition A2 holds. 

(1) In the scalar case {N = 1), De Giorgi-Moser-Nash theory implies A2 with 
fiQ = fj-oin, A), i?c = <», and Bq = i?o(«, A). 

(2) If n = 2 and the coefficients are time-independent, then A2 holds with 
A'o = A^o(A), Rc = GO, and Bq = Bo{N,X). It is a consequence of the well 
known theorem of Morrey [531 PP- 143-148]; see [HI Theorem 3.3] for the 
derivation of A2 from Morrey's estimate. 

(3) Let A"^ be scalar functions on R"+^ such that there is a constant Aq > 
such that 

for arbitrary column vectors t] e M" with elements ^a, Va- We denote 
^= sup ^ 2 L^^(A)-^«^(A)<5, 



,i,j = l a,/3 = l 



Then, there exists S'q = S'o{n,Xo) such that if < <foi then A2 is satisfied 
with ^0 = ^J'o{n, Ao), Rc = co, and Bq = Bo{n, N, Aq); see ^ Lemma 2.2]. 

(4) For a measurable function / defined on R"+^, we set for p > 
^^p = ^^p{f)= sup sup-/- f \f{y,s) - fx,r{s)\dyds 

where fx.r{s) = f f{y,s)dy. 

J B^(x) 

Wc say that / belongs to VMO^: if limp^o'^p(/) = 0. If the coefficients 
belong to VMO^^ class (or their BMOa; modulus Wp arc sufficiently small) 
then A2 holds with Rc, and yUoi Bq depending on n, N, A, fl and the BMOx 
modulus ujp of the coefficients. If is bounded, then we can take Rc = 
diamfi; see [HI Lemma 2.3]. 



4.1.3 Examples for A3 

Finally, we give examples when the condition A3 holds. 

(1) Let = 1 and fl be an extension domain for which there exists 9 > 
and Rm e (0,oo] such that (|3.20p holds for all x e and < i? < Rm- 
Then, A3 holds with Rm and Bi = Bi{n, A, Jl, 9). This follows from the De 
Giorgi-Moser-Nash theorem. 

(2) If n = 2, the coefficients are time-independent, and is a bounded Lipschitz 
domain or a special Lipschitz domain, then A3 holds with Rm = diamJl 
and Bi = Bi{N, X,n). The proof is essentially the same as that of [l3l 
Lemma 4.4]. 
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(3) If the coefficients belong to VMO^; class (or with sufficiently small BMOa; 
modulus ujp) and Q, is bounded domain (or Lipschitz domain with small 
Lipschitz constant), then A3 holds with Rm = diamf2 and Bi depending on 
n, TV, A, fi, and BMO^: modulus ixjp of the coefficients; see [TUl Lemma 6.1] 
and m Corollary 4.11]. 



4.2 Elliptic Neumann function 

Let us consider elliptic differential operator of the form 



Lu = - 

ox, 



-f^"^(x)^), 



where A°'^{x) are N x N matrices whose elements a^^ (x) are bounded measur- 
able functions satisfying (|1.2p . and its adjoint operator L* defined by 



L*u = -^(A^"ixfp-] 
dxa \ dxpj 



Let Vl cz be an extension domain such that |ri| < co. We consider Neumann 
boundary value problem 

fLu = f in f2, 
du/di^ = on dil. ^^'^^ 

Given / e (VF^'^(r2)*)^, we say that u is a weak solution in W^^'^^ft)^ of the 
problem (|4.2p if we have u e W^^'^(fi)^ and 

By the standard elliptic theory, we have a unique weak solution in W^''^{H)^ 
of the problem (|4.2p provided / satisfies the compatibility condition {/, 1) = 0. 
We say that an N x N matrix valued function G{x, y) is the Neumann function 
of L in 51 if it satisfies the following properties: 

i) G{-,y) 6 Wl^^in)'^^ and G{-,y) e W^-'^{n\Br{y))'^^ for all y e n and 
r > 0. Moreover, J^^ G{x, y) dx = 0. 

ii) LG{-,y) = {Sy - in fl, dG{-,y)/dv = Q on dn for all y e SI in the 
sense that for 4> 6 Cf{il)^ satisfying J^4> = 0, we have 

l^^a'^fDpG,,{;y)D^cP^ = cp''{y). 

iii) For any / e C^-'{d)^ satisfying Jq f = 0, the function u defined by 

u{x) = [ G{y,xff{y)dy (4.3) 
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is the weak solution in W^'^(f2)^ of the problem 

L*M = / in rj, 

4.4 

dujdv* = on (5f]. ' 

Note the property iii) and the requirement G{x,y) dx = give the unique- 
ness of the Neumann function. We warn the reader that the above definition is 
slightly different from the one given in a recent article [10] , where it is assumed 
that n ^ 3 and d^l is Lipschitz so that the boundary trace of a W^''^{rt) function 
is well defined. In }10| , the Neumann function satisfies the normalization condi- 
tion J\,jj G(a;, y) dSx = instead of G{x, y) dx = 0. Wc point out that they 
differ by a function U, which is a weak solution in W^'^{n)^ of the problem 

LU = in n, 

dU/dv=—j^I on dVl. 

We treat only the two dimensional case below because the case when n ^ i 
is studied in [10] extensively by a different but more direct method. We also 
mention that in a recent article |30| . the Neumann function is constructed in 
BMO(ri), where i7 is a bounded Lipschitz domain in R^. 

Theorem 4.5. Let 51 c be an extension domain such that < co. Then, 
there exists the Neumann function G{x, y) of L in H. It is continuous in {(x, y) e 
fl X fl: X ¥= y} and for x, y satisfying < |a; — y| < i?/2, where R = d^ /\ dy, we 
have 

\Gix,y)\^C{l + g^R-' + HR/\x-y\)), (4.6) 
where C = C{N,X) and g is the same constant as in (|5.42l) . Also, we have 

G{y,x) = G*{x,yf (4.7) 

where G*{x,y) is the Neumann function of the adjoint L* in f2. Moreover, if 
Q is a Lipschitz domain, then for all x,y e Q, we have 

\G{x,y)\tiC{l + \ii{d/\x-y\)), d = diamQ, (4.8) 

where C depends on the Lipschitz character of dQ, dianiH as well as N and A. 

5 Proofs of main Theorems 

5.1 Proof of Theorem 13. 9t Case when < co. 

In the proof, we denote by C a constant depending on n, N, X, no, Bq unless 
otherwise stated; if it depends also on 7, it will be written as C~^, etc. We 
begin with constructing mollified Neumann Green's functions. Fix a function 
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$ G C^{Bi(0)) satisfying < $ 2 and 4„ 4> = 1. Let Y = {y,s)eQ be fixed 
but arbitrary. For e > 0, we define 

$,,,(x) = e-"$((.T-y)/e). 

For any T > s, let = v^-y.k be a unique weak solution in ';^^'°(f2 x (s,T))^ 
of tlie problem 

^u = in nx (s,r) 
du/di' = on dn X (s,r) (5.1) 
u = ^y i^Bk on X {s}, 

where is the fc-th unit vector in R^; with an aid of the inequality p.Sp . which 
was a consequence of the condition Al, the unique solvability of the problem 
(|5.ip follows from the Galerkin method described in [551 §111.5]. By setting 
Vt{x,t) = for t < s and letting T — > cx), we may assume the is defined on 
the entire Q. Then, by the energy inequality, we obtain 

\\\v4\Q^C\\%4LHn)^Ce-"/'- (5.2) 

Thus t)e belongs to "^(Q)^. We define the mollified Neumann Green's function 
N\;Y) = (7V|,(-,r))^^^^^ for ^ by setting 

Next, for a given / = (/i, . . . , f^)^ e "€^{^1 x (a, 6))^, where a < s < 6, let it 
be a unique weak solution in "^^'"(il x (a, 6))^ of the backward problem 

= f in ilx (a, b) 
du/di^* =0 on 5r2 X (a, 6) (5.3) 
u = on fl X {b}. 

By setting u(x, t) = for t > 6 and letting a — oo, we may assume that u is 
defined on the entire Q. Then, it is easy to see that we have the identity 

/ %,,{x)u''ix, s) dx = f[ Ntk{X, Y)f\X) dX. (5.4) 
Jn JJ Q 

If we assume that / is supported in Q = (5j(Xo), where < R < dy, then the 
energy inequality and Holder's inequality yields 

II|w|IIq «^ C'll/ll^2(„ + 2)/(„ + 4)(Q)ll"ll^2(„ + 2)/„(Q)- 

Then, by the inequality p.Sp applied Q = Q'^{Xo) = Br{xo) x {to, to + i?^), we 
get 

II|w|||q<C||/||^,,„^,,,^„,,,(q). (5.5) 

By utilizing (|5.5p and the condition A2, and proceeding as in |5J Section 3.2], 
we find that u is continuous in jQ = Q+y^(Xo) and satisfies (see [SI Eq.(3.15)]) 

Ho;iQ<C^fl'll/ll^.(Q)- (5-6) 
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If B,{y) X {s} c ig, then by (gH) and dSTB]) . we obtain 



0;4( 



Therefore, by duahty, it foUows that we have 

\\N'i-,Y)\\^^^Q^^CR', (5.7) 

provided < R < dy and B^{y) x {s} c = (5+^^(Xo). For X e Q such that 
< d:= \X - Y\g^ < dy/Q, if we set r = d/3, Xq = (y, s - 2^^), and R = 6d, 
then it is easy to see that for e < d/3, we have 



and also that = v^.^y^k satisfies = in (X). Then, we derive from 

the condition A2 that (see [8j Lemma 2.6]) 



By combining (|5.7p and (|5.8p . we obtain that 

|t;,(X)|<Cr-«-2|I^,||^^^^_(^^, 

That is, for any X,Y e Q satisfying < \X — < dy/6, we get 

\N'{X,Y)\iiC\X -Y\-^, \fe< l\X -Y\.g,. 
Next, notice that we have 

/ N'{x,t,y,s)dx = H{t-s)I; H{t) = l[o,a,){t). 
Jn 



We define 



so that 



N'{x, t, y, s) = N'{x, t, y, s) - —Hit - s)I 



N%x,t,y,s) dx = 0. 
Therefore, by (|3.7p . for any — G0sga<5sgcx), we have 

l|A^'(-,y)ll^.(„ + .,/..(f^xK'')) < 7ll|A^'(-,i^)ll|ox(a,6). 

Lemma 5.13. For any Y e Q, < R < dy, and e > 0, we have 

\\\N^{;Y)\\\Q^Q^^y)^CR-/', 



(5.8) 



(5.9) 



(5.10) 
(5.11) 

(5.12) 

(5.14) 
(5.15) 
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Also, for any Y e Q and e > 0, we have 

\{X e Q: \N'{X,Y)\> t}\ iiC^T-'^, \fT>dy'', (5.16) 
\{X e Q: \D^N'{X,Y)\> t}\ iiCT-'^, Vr > (5.17) 

Furthermore, for any Y e Q, < R < dy , and e > 0, we have 

\\N%;Y)\\^^^Q^^y))^C,,,R-^ + ^-+^y^, Vp6[l,i^), (5.18) 
PAr^(.,y)||^^(Q,(y)) Cpi?— vpe [l,^). (5.19) 

Proof. By JOl), (EH), and the energy inequality, we get (see [HI Eq. (3.20)]) 

\Wi-^y)lQ\Qn(Y) < CR-^^I\ 0<R< dy/G. 

from which and (|5.14p follows readily since 1 57 1 5= Ci?" and dy is comparable to 
Jy/6. In fact, if we fix a C 6 'rf^ {Q niX)) , where < i? < dy/d, such that 

0«;C<1, C= 1 onQ^/2(y), \D,Cf + \Ct\ ^ 16R-^ (5.20) 
then we have (see [HI Eq. (3.22)]) 

111(1 -c)A^^(•,>^)lllQ<ci?-"/^ 

and thus by (|5.10p and the properties of (, we get 

111(1 - ON'{;Y)\\\q < C(i?-"/2 ^ i^n/2|^|-l + < Ci?-"/^^ (5.21) 

To get (I5T^ . set Q(i) := f7 x (s + i?^ (») and Q(2) := {n\Bji{y)) x (s, s + i?^) 
and note that by (|5.12p and (|5.14p we have 

ll^^(^>^)ll^.™,.(Qa,) < C'7|||iV^(-,r)|||g^,, ^ C7i?^"/2 (5.22) 
and similarly, by p.Sp and (|5.2ip . we have 

>^)ll^.(„+.,/„(Q(.,) < (27 + 1)^ 111(1 - C)A^^(-, y)lllnx(.,.+fl^) 

=s; C(27+ l)^i?-"/2. (5.23) 

Therefore, by combining ()5.22p and (|5.23p . we obtain (|5.15|) . We derive (|5.18p 
and (|5.19p . respectively, from (|5.16p and (|5.17p . which in turn follows respec- 
tively from (|5.15|) and 1)5.14^ : see [H Lemmas 3.3 and 3.4]. ■ 

Lemma 5.24. Let {ukYf^j^ be a sequence in "^(Q). // supj, |||ufc|||g ^ A, then 
there exist a subsequence {ukjj'jLi — {^^IfcLi '^"'^ ^ ^(Q) with |||Mj||Q ^ A 
such that Ukj u weakly in W2'^{^ x (a, b)) for all — oo < a < b < (X>. 

Proof. See Lemma A.l]. ■ 
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The above two lemmas contain all the ingredients for the construction of a 
fmiction N(-,Y) such that for a sequence tending to zero, we have 

7V^-(.,r) ^ N{;Y) weakly in <'°(grfv(r))^', 
(l-C)iV'^(-,r) — (l-C)iV(-,r) weakly in #2^'°(17 x (-T,r))^', VT > 0, 

where 1 < q < and ( is as defined (|5.20p with R = dy/i. It is routine 
to verify that N{-,Y) satisfies the same estimates as in Lemma [5.131 see |51 
Sec. 4.2]. Note that ((glT|) implies N{x, t, y, s) dx = 0. We define N{-,Y) by 

Nix, t, y, s) = N{x, t, y, s) + ^^H{t - s)I. (5.25) 

Then, by the previous remark that it is clear that N{X, Y) satisfies the local 
estimates in Remark 13.161 and thus the property a) for the Neumann Green's 
function defined in Sec. 12.41 We shall now show that N{X, Y) also satisfies the 
properties b) and c) so that N{X, Y) is indeed the Neumann Green's function. 
To verify the property b), let us assume <p = is supported in 

n X {-T,T), where -T < s <T, and note that 

N'-{;Y) ^ N{;Y) weakly in W^'\Qi^.{Y)f\ (5.26) 
{I- QN^^^i-.Y) ^ {I- C)N{-,Y) weakly in #'2^^°(17 X (-T,r))^'. (5.27) 

Since fcth column of N^{-,Y) is the weak solution in x (s,r))^ of the 

problem (|5.1[) . we find that 

= fj^ [-Nli:{X,YW,{X) + atfDpN;^^{X,Y)D^cl>\X)] dX 

By writing (p = rjcj) + {1 — 'q)4>, where 77 e ^^{QiIyO^)) satisfying 77 = 1 on 
QdY 12(^)1 using ()5.26|) . (|5.27p . and taking ^ ^ cx) in the above, we get the 
identity (|2.4p : sec [SI p. 1662] for the details. To verify the property c), let us 
assume that / is supported in f2 x (a, b), where a < s < b and u be the unique 
weak solution in 1^^'°(ri x (a,b))^ of the problem (|5.3p . By setting u{x,t) = 
for t > b and letting a — > —00, we may assume that u is defined on the entire 
Q. Then, similar to ()5.4p . we have 

[ ^y.,^Xx)u\x,s)dx= [[ Nl-{X,Y)f\X)dX. 

By the condition A2, it follows that it is locally Holder continuous in Q; see 
the remark we made in deriving (|5.6p . By writing / = C/ + (1^C)/ ^^nd using 
(|5.26|1 . ()5.27p . and taking the limit /.j ^ <», we then get 

u{Y)= [ N{X,Yff{X)dX. 
JQ 
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Therefore, we have u = u and thus the property c) is verified. 

It is clear from the construction that N{x,t,y, s) = if t < s. For any 
X = {x, t) e Q, we define the moUified Neumann Green's function of 
with a pole at X by letting its l-th column to be the unique weak solution in 
1^^'°(r2 X (T,t))^, where T < t is fixed but arbitrary, of the backward problem 

if *i6 = in 17 X (T, t) 
dujcv* = on Pf] X (T, t) 
u = ^x,eei on X {t}. 

As before, we extend it to the entire Q by letting it to vanish on ^l x (t, co). By a 
similar argument as above, we obtain the Neumann Green's function N*{-, X) 
of .^f* that satisfies the natural counterparts of the properties of the Neumann 
Green's function for Note that the condition A2 together with the esti- 
mates i), ii) listed in Remark l3.16l and its counterparts imply that N{-,Y) and 
N*{-,X) are locally Holder continuous in Q\{Y} and Q\{X}, respectively. Us- 
ing the continuity discussed above and proceeding similar to [51 Lemma 3.5], we 
find that 

N{Y,X) = N*{X,Yf, \/X,YeQ, X ^ Y. (5.28) 
Moreover, similar to |S1 Eq. (3.44) and (3.45)], for t > s, we have 

N^ix^t^y^s) = / ^y^,{z)N{x,t,z,s)dz, 
Jn 

which justifies why we call it "mollified" , and 

lim N'{x,t,y,s) = N{x,t,y,s). (5.29) 

By (|5.28p and the counterpart of the property c) in Sec. 12.41 we sec that u 
defined by the formula ([QUI) is a weak solution in 'r^'°{Q)^ of (|XTT|) . 

We now prove the identity (|3.12p for the weak solution in "^^'"^(il x (s, co))^ 
of the problem (|3.13p . Similar to (|5.4p . we have 

<Px,,iy)u''{y,t)dy = / Ntk{y,s,x,t)iP\y)dy, 
Jn 

where iV^ (• , X) is the mollified Neumann Green's function of ^* . The condition 
A2 implies that u is continuous in x (s, cx)), and thus we have 

lim/ '^.x^,{y)u^{y,t)dy = u''{x,t). 

On the other hand, by (|5.29p and the counterparts of (|5.14p . together with the 
dominated convergence theorem, we get 

lim / N^,{y,s,x,t)r{y)dy= f NUy,s,x,t)ij\y) dy. 
^^°Jn Jn 

Then, the identity ((3?T2)) follows from ((5?28| . Finally, we obtain ([3TT4)) similar 
to ([5^ and get (P?T5)) from ([XTi)) and the condition A2. ■ 
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5.2 Proof of Theorem 13. 9t Case when = co. 

The proof for the case when \U\ = co is almost identical to the case when 
|n| < 00. As a matter of fact, it is even more simple since we do not need to 
introduce N'^{X,Y). From the condition Al, we obtain 

ll^'(-'^)ll^2(„ + 2,/„(Ox(a.6)) < 7\\\N'i-,Y)\\\n>,ia.b). 

By using the above inequality instead of (j5.12p and follow the same augment as 
in Section [5.11 we get the same conclusion. ■ 



5.3 Proof of Theorem IST^ 

We adopt the argument in [TH] , which is in turn based on [TU [TS] . Let -0 be a 
bounded Lipschitz function on M" satisfying \DiIj\ Al a.e. for some M > be 
chosen later. For i > s, we define an operator Pf_^f on L^{n)^ as follows: For 
/ e L2(j7)Af and T > t, let u be the unique weak solution in -f^^^^Q x (s,T))^ 
of the problem 

\^u = in fix {s,T), 
du/df = on oil X (s,r), (5.30) 
u = e^'f'f on fl X {s} 

and define Pf^tfix) ■= e'^(^)M(a;, t). Then, by (l3T2l) . we find 

Pf-.tf{x) = e'/'(-) / iV(x, t, y, s)e-^'^y^f{y) dy. (5.31) 
Jn 

Then, as in [SI Sec. 5.1], we derive 

\\Pttf\\LHn)^e''''^'-^^\\fh.^n), (5.32) 
where i? = A~^. We set p = -s/t — s a Rm and use the condition A3 to estimate 

e-'^<^^^\Pttf{^)\''=H^.t)\' 

,2 -(n+2) f f I ,/ , x|2, 



^Btp-^-+'> / \u{y,T)\Uydr 



^52^-(n+2) f f e-^^^y)\Pf_^,m\^dydr. 
Thus, by using (|5.32p . wc derive 

\Pttf(x)\' ^ BIp—' f f e^^^^^-^^(y^\Pt^f{y)\^dydr 

^Blp—^f f e^'>'^\P:trfiy)\'dydr 

Jt-p^ JnnBp{x) 
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We have thus obtained the following estimate for Pf^t'- 

WPLjWu^n) < i3ip-"/'e''™^^('-^)||/|U.(o). (5.33) 

We also define the operator Qf^^ on L^(ri)^ for s < t by setting Qf^sg{y) = 
e~'^^y'>v{y, s), where v is the weak solution in "/^^'"(fi x (T,t))^ , where T < s, 
of the backward problem 

= in n X {T,t), 
du/dv* = on dnx {T,t), (5.34) 
u = e^^g on X {t}. 

Similar to (|5.33p . we have 

WQtsah^Hn) < Sip-"/^eP^-'+"''^'*-^^ll9llL^(n). (5.35) 
Note that by (|O0)) and ([OH) we have 

/ {Ptj)-g = [ f-{Qts9)- 
Jn Jn 

Therefore, by duality, (|5.35p implies that for all / e C^(r2)^, we have 

WPttfhHn) ^ SiP^"'''e''^^+^^^^(*-)||/|U.(o). (5.36) 
Now, set r = (s + t)/2 and observe that by the uniqueness, we have 

Then, by noting that t — r = r — s = {t — s)/2 and p/\/2 -sjt — r a Ru p, 
we obtain from ([031 and ([061 that for all / e (f^)^, we have 

where C = 2^I^B\. For all 6 with x ¥= y, the above estimate combined 
with (|5.3ip yields, by duality, that 

e'A(^)-'A(!/)|7v(x,i,2/,s)| ^ Cp-" 62"^^+''^^' (5.37) 
Let ip{z) ■■= MtpQ^z — y\), where i/'q is defined on [0, go) by 



r ii r ^ \x — y\ 

\x - y\ if r>\x- y\. 



V'o(f) 

Then, -0 is a bounded Lipschitz function satisfying \Dip\ M a.e. Take M = 
\x — y\/2i!)(t—s) and set r = \x — y\/^yt — s. By (|5.37|) and the obvious inequality 
pl\ft^^ ^ 1, we have 

|Ar(x,t,y,s)| Cp-" exp(r/i5 - rV4i?). 
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Let A = A{-d) = A{X) be chosen so that 

exp(r/i9 - T^/M) s: ^cxp(-rV8i?), Vr e [0, oo). 
If we set K = 1/8?? = then we obtain 

\N{x, t, y, s)\ < Cp'" exp {-k\x - y\^/{t - s)] . 
where C = C(n, N, \,Bi). We have thus proved FinaUy, ((3?22l) gives 

I Ar(x, t, y, s)K C {1 V ((t - {t - exp - y\^/{t - s)} , 

which clearly implies ifO<<-s<r. ■ 

5.4 Proof of Theorem IST^ 

Notice that (|3.25|) implies that 

\N{X,Y)\^C\X-Y\-^' iiQ<\t-s\<Rli, (5.38) 

where C = C{n, N, A, B2,k). Then, by the energy inequality (see [HI Eq. (3.21)]) 
and noting that Rm/'2 and Rm are comparable to each other when Rm < oo, 
we derive from (|5.38p that for < ?' < Rm, we have 

lll^(-,nills\Q.(y)<C'r-"/2. (5.39) 

Let X e Q and < i? < be given. Without loss of generality, we assume 
X = Q and write Qn = Qr{0), etc. Let it be a weak solution of 

= in Q = Q+ n Q, du/di''' = Q on S = Q% n dpQ. 

We will show that u is locally bounded in ^Q, where we use notation aQ = 
QaR 2' ^'^^ satisfies the estimate (|3.18p . The other case can be treated in 
the same fashion. Let w = rju, where rj e '^^''{QRft) is such that 

77 ^ 1, 77 = 1 on ggfl/s, \D,if + |77t| 32/i?^. 

Then, w becomes a weak solution of 

^*-w = -Tjtu - {A^°')'^DpuDaii - Da {{A°'^)^DfjT]u) infix (0, <»), 
dw/dv* = {AP°')'^ Dprjuva on dVt x (0,cx)). 

For Y = {y,s) e jQ, let N'^{-,Y) be the mollified Neumann Green's function 
as constructed in the proof of Theorem 13.91 Then similar to (|5.4p , we have 



$y,e(a;)77(a;,s)u''(x,s)dx = - / 77t7i'iV4(-, F) 

a^tDpuW^r^N^^{;Y)+ [ a^^^ D pi^uW aNt^{- ,Y) . 
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Since \X — Y\,^ > R/8 for X e Q r\ supp{\Dj;r]\ + |»7t|), the integrals on the right 
hand side are in fact integrals over ^Q\Q r/sO^) ■ Then, we take limit e to zero 
in the above by using (|5.26p . (|5.27|) to get 

^HY) = - [ Vtu'N,k{;Y) - [ affN,k{-,Y)D^u'Dpr, 
HQ HQ 

+ f afDpNjk{-,Y)u'D^r^ =: h+h + h- (5.40) 
HQ 

Denote Au{Y) :=Qn {Q3R/4{Y)\Qji/s{Y)) zd ^Q\QR/siY). By the properties 
of 77, Holder's inequality, and (|5.39p . we estimate 

I/3I < Ci?-ip,Ar(,r)||^,(^,(y)) ||w||^,(iQ) < 

Also, by using (|5.38p and the energy inequality, we estimate 

CR-'\\N{;Y)\\j^^^^^(^Y)) ll^«ll^.(iQ) < Ci?-("+2)/2||M||^,(Q). 

By combining above estimates for /i , /2 , and I3 , we conclude from (|5.40p that 

where C = C{n, N, A, B2, k). Then, we obtain the condition A3 from the above 
inequality by a standard covering argument. I 

5.5 Proof of Theorem 14.51 

The idea of proof is the same as that of [T31 Theorem 2.12] but wc reproduce 
some key steps of the proof for the completeness. Let K{x, y, t) = N{x, t, y, 0), 
where N{x, t, y, s) is as in (|5.25p and N{x, t, y, s) is the Neumann Green's func- 
tion oi J^f = 8/ 8t + L in Q. We define the Neumann function G{x, y) by 

G{x,y) := / K{x,y,t)dt. (5.41) 

Since f2 is an extension domain and |f2| < <», there is a constant g = g{^l) such 
that for all It e W^^'^(r2), we have 

ll"llL^(n) < g\\Du\\L2(^n)- (5.42) 

Then by [131 Lemma 3.12], wc find that the integral in (|5.4ip is absolutely 
convergent for x ^ y and thus G{x, y) is well defined. Also, it is clear that 
J^G{x,y) dx = 0. We define the Neumann function G*{x,y) of L* similarly. 
Then, by wc find that (see [H Eq. (3.21)]) 

G*{x,y)= N*{x,-t,y,0)dt= K{y,x,tfdt=G{y,xf, 



23 



where N*{x,t,y,s) = N*{x,t,y,s) — -^H{t — s)I and N* {x,t,y, s) is the 
Neumann Green's function for ^* = —d/dt+ L* . Let us denote 



K{x,y,t) = / K{x,y,s)ds 
Jo 



Lemma 5.43. The following estimates hold uniformly for t > and y e fi. 

a) \\K{;y,t)\\L.(Biv,d,)) *S Cp{dl + g^)dl^P-\ Vpe [1,2). 

b) \\D,K{;y,t)\\L,^B(y^,^)) < Cp{dl + Q^) dy , Vp £ [1,4/3). 

c) \\k{;y,t)\\LHn\B(y,r)) ^ C{r^ + g^)r'^/\ \/re{0,dy]. 

d) \\D^k{;yML-(n\B(y,r)) ^G{r^ + g')r-\ Vre(0,d^]. 

Proof. See [HI Lemma 3.23]. ■ 

Wc have to show that G{x,y) defined by the formula (|5.4ip satisfies the 
properties i) - iii) in Section 14.21 We focus on iii) , which may seem less clear 
than i) or ii). For any / e C"(r2) satisfying J^^ f{x) dx = 0, let u be defined by 
Note that 

v{x,t) := [ k{y,x,tff{y)dy = [ [ K{y,x,sff{y)dyds 
Jn JQJn 

is absolutely convergent by Lemma 15.431 and 

limv{x,t)= / G{y,xff{y)dy = u{x), (5.44) 

vt{x,t)= [ K{y,x,tff{y)= [ N*{x,-t,y,0)f{y)dy, 
Jn Jn 

where we have used the assumption / = in the last step. Therefore, for 
any T > 0, w ■■= Vf is the weak solution in 'f^^'^i^ ^ (0, T))^ of the problem 

Wt + L*w = in r2 X (o,r) 
dw/dv* = on dn X (o,r) 
w{-,0) = f on n, 

Then, as in the proof of [131 Lemma 3.12], for t > 0, we have 

\\vt{;t)U2^n)^Ce-^^'"'\\fU2^n), (5.45) 

which implies Vt{-, t) e L^{il,)^ . Also, by using N*{x, t, y, s) = N*{x, t—s, y, 0), 
we find 



v{x,t) = / N*{x,-t,y,s)f{y)dyds, 
J -t Jn 
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and thus, for T > 0, i; is the weak solution in "/^^'°((0, T) x $7)^ of the problem 

' vt + L*v = f in nx (0,T) 

dv/dv* = on dflx (0,r) (5.46) 
v{-,0) = on n. 

Then, by using ([QS]) and ^M^, we obtain (see [H Eq. (3.47)]) 

\\DM;t)\\LHn)^C\\f\\L2^n), a.c. i > 0. 
Therefore, by (|5.42p for a.e. < > 0, we have 

\\v{-,t)\\w^-2(n) ^ ^11/11^2(0). 

Then, by the weak compactness of the space W^'^{^1), we find that there is 
a sequence tm oo such that v{-,tm) u weakly in VF^'^(ri)^ for some 
u 6 W^'^{n)'^. By ([5^ . we must have u = u and similar to [HI Eq. (3.48)] 
and the argument after it, we find m is a weak solution of the problem (|4.4p . By 
using (|4.7p . it is easy to see that J^u = 0. We have thus verified the property 
iii). By repeating essentially the same proof of [SI Theorem 2.12], we find that 
G(.T, y) satisfies the properties i) and ii) as well, and that it has the logarithmic 
bound ((46| . 

If is a Lipschitz domain, then by modifying [T31 Lemma 4.4], one can 
show that ^ = d/dt + L and ^* = —d/dt + L* satisfy the condition A3 with 
Rm = diam 51 = d and Bi depending on A, N and Lipschitz character of (557; 
see Section l¥. 1.31 By (|5.38p . for X = {x,t) & Q satisfying \t\ < (P , we have 

\K{x,y,t)\^C\X-Y\-^, Y = {y,Q), (5.47) 

and by (|5.39p . for < r < c?. we have 

lK{-,y)lQ\Q(Y,r)^Cr~\ 

Similar to [HI Eq. (3.59)], for < r < d and t^2r'^, we have 

\K{x,y,t)\ < Cr-'^e~^e-'{t-2r')^ (5 ^g-. 

We set r := i min(gi, d). U < \x — y\ ^ r, then by (|5.4ip . we have 

r\x-y\^ r2r^ /.GO 

|G(x,y)l^/ +/ +/ \K{x,y,t)\dt=:h+h + h. 

Jo J\x-y\^ J2r^ 

It then follows from (|CT7| and (fOSj) that 

|2 

f\x-v\ 

/i C / \x-y\-^dt ^ C, 
Jo 

/■2r^ 
r^dtiiC + C\n{r/\x-y\), 

J2r^ 
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Combining all together we conclude that ifO<|x — y|^r=| min(gi, d), then 
\G{x,y)\ < C (1 + (e/r)2 + \n{r/d) + \n{d/\x - y\)) . (5.49) 
li \x — y\ ^ r = ^ niin(£i, d), then by (|5.47p and (|5.48p we have 

|G(x,y)| < / + / \Kix,y,t)\dt 

Jo J2r^ 

r-^ + C r-^e-^^'"^'-^'"Ut^C+CQ'r-^. (5.50) 

By combining (|5.49p and (|5.50p . and using d < cx), we get (|4.8p . ■ 
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